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Abstract 

The J/^> mass shift in cold nuclear matter is computed using an effective Lagrangian approach. 
The mass shift is computed by evaluating D and D* meson loop contributions to the J/^f self-energy 
employing medium-modified meson masses. The modification of the D and D* masses in nuclear 



matter is obtained using the quark-meson coupling model. The loop integrals are regularized with 
dipole form factors and the sensitivity of the results to the values of form-factor cutoff masses is 



> 

O 

p,i investigated. The J/\P mass shift arising from the modification of the D and D* loops at normal 



nuclear matter density is found to range from —16 MeV to —24 MeV under a wide variation of 
values of the cutoff masses. Experimental perspectives for the formation of a bound state of J/^l 
to a nucleus are investigated. 
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I. INTRODUCTION 

A new era of nuclear matter research is envisaged with the 12 GeV upgrade of the 
CEBAF accelerator at the Jefferson Lab in the USA and with the construction of the FAIR 
facility in Germany. These new facilities will have the exciting potential of implanting 
low-momentum charmonia and charmed hadrons in an atomic nucleus, like the J/^f and i/j 
mesons and heavy-light charmed mesons such as D and D*. While at JLab charmed hadrons 
will be produced by scattering electrons off nuclei, at FAIR they will be produced by the 
annihilation of antiprotons on nuclei. There are several reasons for the excitement, one of the 
main ones being the opportunity of studying the poorly understood low-energy excitations 
of gluon degrees of freedom. An example where these excitations play an important role is 
the propagation of charmonia in matter. Since a charmonium state does not have quarks in 
common with the nuclear medium, its interactions with the medium necessarily involve the 
intervention of gluons. Basic interaction mechanisms discussed in the literature have been 
the excitation of QCD van der Waals forces arising from the exchange of two or more gluons 



between color-singlet states 



2] , and the excitation of charmed hadronic intermediate states 



with light quarks created from the vacuum [3|, |4j • 

Another interesting challenge is to study the properties of charmed D and D* mesons in 
medium. The chiral properties of the light quarks that compose these mesons are much more 
sensitive to the nuclear medium than their companion, heavier charm quarks and therefore 
they offer the unique opportunity of studying phenomena like the partial restoration of chi- 
ral symmetry in nuclear matter. Motivated by such considerations, some very interesting 
phenomena involving these mesons have been predicted. Amongst these we mention the 
possible formation of D(D) meson-nuclear bound states [5j, enhanced dissociation of J/^f 
meson in nuclear matter (heavy nuclei) [6J], and enhancement of the D and D meson produc- 
tion in antiproton-nucleus collisions [7J. Ref. [8j presents a recent review of the properties 
of charmonium states and compiles a fairly complete list of references on theoretical studies 
concerning a great variety of physics issues related to these states. On the experimental 
side, one of the major challenges is to find appropriate kinematical conditions to produce 
these hadrons essentially at rest, or with small momentum relative to the nucleus, as effects 
of the nuclear medium are driven by low energy interactions. 

The original suggestion [2| was that QCD van der Waals forces arising from multiple gluon 



exchange would be capable of binding a charmonium state by as much as 400 MeV in an 
A = 9 nucleus. The estimate was based on a variational calculation using a phenomenological 
ansatz for the charmonium-nucleus potential in the form of a Yukawa potential. Along the 
same lines but taking into account the distribution of nucleons in the nucleus by folding the 
charmonium-nucleon Yukawa potential with the nuclear density distribution, Ref. 9| found 
a maximum of 30 MeV bin ding energy in a large nucleus. A somewhat more QCD-oriented 
estimate was made in Ref. lOj. Using a lowest-order multipole expansion for the coupling 
of multiple gluons to a small-size charmonium bound state jl], it is possible to show on 
the basis of the operator product expansion that the mass shift of charmonium in nuclear 
matter is given, in the limit of infinitely heavy charm quark mass, by an expression similar to 
the usual second-order Stark effect in atomic physics, which depends on the chromo-electric 
polarizability of the nucleon. Using an estimate [1] for the value of this polarizability, the 
authors of Ref. [lOj obtained a 10 MeV binding for J/ty in nuclear matter. On the other 
hand, for the excited charmonium states, a much larger binding energy was obtained, e.g. 
700 MeV for the excited charmonium state tp'(2S), an admittedly untrustworthy number. 
Following this same procedure, but keeping the charm quark mass finite and using realistic 
charmonium bound-state wave-functions, Ref. |4j found 8 MeV binding energy for J/ty in 
nuclear matter, but still over 100 MeV binding for the charmonium excited states. While 
an increase in the QCD Stark effect is expected for excited states (because of their larger 
size), the extreme values for the binding energies for these states found in the literature are 
widely considered to be unrealistic. The source for such an overestimate is attributed to the 
breakdown of the multipole expansion for the larger-sized charmonium states. 

There are some other studies on charmonium interactions with ordinary hadrons and 



nuclear matter, in particular involving the J/^> meson. QCD sum ru 



es studies estimated a 



J/ty mass decrease in nuclear matter ranging from 4 to 7 MeV [llMl3|. while an estimate 
based on color polarizability [l4J gave larger than 21 MeV. In addition, there are studies 
of the charmonium-nucleon interaction and of J/\l/ dissociation cross sections based on a 
one-boson exchange model (l5j, effective Lagrangians [16l Il7| and the quark-model [18]. In 
Ref. [l9j the charmonium-hadron interaction was studied in lattice QCD. 

A first estimate for the mass shifts of charmonium states (we denote charmonium states 
generically by ip) in nuclear medium arising from the excitation of a pair of D and D* mesons 
- see Fig. [I] - was performed in Ref. [4 1 . Employing a gauged effective Lagrangian for the 



coupling of D mesons to the charmonia, the mass shifts were found to be positive for J/\l/ 
and ip(3770), and negative for -0(3660) at normal nuclear matter density po- These results 
were obtained for density- dependent D and D masses that decrease linearly with density, 
such that at po they are shifted by 50 MeV. The loop integral in the self-energy (Fig. [T]) 
is divergent and was regularized using form-factors derived from the 3 Po decay model with 
quark-model wave functions for ip and D. The positive mass shift is at first sight puzzling, 
since even with a 50 MeV reduction of the D masses, the intermediate state is still above 
threshold for the decay of J/ty into a DD pair and so a second-order contribution shoulc 
be negative. As we shall explain below, this was not realized in the calculation of Ref. |4] 
because of the interplay of the form factor used and the gauged nature of the interaction. 
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FIG. 1. DD-loop contribution to the J/\& self-energy. We include also DD* and D*D* contribu- 
tions. 

In the present paper we reanalyze the mass shift of J/^f in terms of the excitation of 
intermediate charmed mesons using effective Lagrangians. In addition to the DD loops, we 
also include DD*, D*D and D*D* loops. The medium dependence of the D and D* masses 



is included by an explicit calculation using the quark-meson coupling (QMC) model 20 ]. 
The QMC is a quark-based model for nuclear structure which has been very successful in 
describing nuclear matter saturation properties and has been used to predict a great variety 
of changes of hadron properties in nuclear medium. A review of the basic ingredients of the 



model and a summary of results and predictions can be found in Ref. 21]. 

The paper is organized as follows. In the next Section we present the effective Lagrangians 
used to calculate the J/^f self-energy and give explicit expressions for the contributions of 
the different intermediate states. In Section Hm we briefly review the QMC description of the 
D and D* mesons in nuclear matter and present numerical results for the density dependence 
of the D and D* masses. A full set of numerical results for the density dependence of the 



J/ty self-energy is presented in Section [TV] We show results for the separate contributions 
of the DD*, D*D and D*D* loops and also investigate the sensitivity of our results to the 
cutoff masses. Our conclusions and perspectives for future work are presented in Section [Vj 

II. EFFECTIVE LAGRANGIANS AND J/* SELF-ENERGY 

We use the following phenomenological Lagrangian densities for the vertices J/^-D and 
J/fy-D* (in the following we denote by ip the field representing J/^)'- 

C fDD = ig^DD V [D {d^D) - (d^D) D], (1) 

/W = 9J£B1 Safi,, (d a ^) \{d»D* v ) D + D (cW)l , (2) 

fn^p l 



+ [(d,r) d: - v faD:)] d^ + d*» [r {d„pi) - (<vn my (3) 

Our convention for the D-meson-field isospin doublets is 

D = (D° D~), D= I " I. (1) 




We note that these Lagrangians are an SU(4) extension of light-flavor chiral-symmetric 
Lagrangians of pseudoscalar and vector mesons. In the light flavor sector, they have been 
motivated by a local gauge symmetry principle, treating vector mesons either as massive 
gauge bosons or as dynamically generated gauge bosons. In the first case, there appear con- 
tact interactions involving two pseudoscalar and two vector mesons. When extended to the 
charm sector, in Eq. (JT]) for instance, there is an additional term of the form 2g"j }DD ijj IJ "ip ^DD . 
In view of the fact that S77(4) flavor symmetry is strongly broken in nature, and in order 
to stay as close as possible to phenomenology, we use experimental values for the charmed 
mesons masses and use the empirically known meson coupling constants. For these reasons 



we choose not to use gauged Lagrangians - a similar attitude was followed in Ref. [22j in a 
study of hadronic scattering of charmed mesons. However, in order to compare results with 
Ref. |4j and assess the impact of a contact term of the form ipipDD, we will also present 
results for the J/^f mass shift including such a term. 



We are interested in the difference of the in- medium, ml, and vacuum, rrif, masses of 



J/*, 



Am = ml — m^, 



(5) 



with the masses obtained from 



ml = K) 2 + £(fc 2 = ml) 



(6) 



Here mB, is the bare mass and S(/c 2 ) is the total J/^S self-energy obtained from the sum 
of the contributions from the DD, DD* and D*D* loops. The in-medium mass, ml, is 
obtained likewise, with the self-energy calculated with medium-modified D and D* meson 
masses. 

We take the averaged, equal masses for the neutral and charged D mesons, i.e. m^o = 
m£>± and m^.o = m£,*±. Averaging over the three polarizations of J/^/, one can write each 
of the loop contributions to the J/\l/ self-energy £/, I = DD, DD*, D*D*, as 



Zi{ml 



3^2 



rfgq 2 F / (q 2 )^(q 2 ), 



(7) 



where -F)(q 2 ) is the product of vertex form-factors (to be discussed later) and the -ft'z(q) for 
each loop contribution are given by 
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4m l/) u; jD 
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U D * 



A(q° = u D * + J7ty) 



wd* - m^/2 u D * + m^/2 



(8) 

(9) 

(10) 



where u D = (q 2 +m|,) 1//2 , uj d * = (q 2 +m|,*) 1//2 , uj d = (uo D +uo D *)/2, £ = for the non-gauged 
Lagrangian of Eq. ([T]) and £ = 1 for the gauged Lagrangian of Ref . 4| , and 



A{q) = J2M<l)i 



(11) 



i=l 



with 



g 2 + (g-^) 2 , [q-{q-k)f 



A(g) = -4g 2 <|4- ^ 
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(12) 
(13) 

(14) 

(15) 



In these last expressions, g and k are four-vectors given by q = (g°, q) and k = (m^, 0). 

III. QUARK-MESON COUPLING MODEL AND D AND D* MESONS IN MAT- 
TER 



In this section we briefly review the QMC description of the D and D* mesons in nuclear 

nn 

matter. Notations and explicit expressions are given in Refs. [51. 1261] . 

The QMC model was created to provide insight into the structure of nuclear matter, 
starting at the quark level [20|, |2l|, |23j. Nucleon internal structure was modeled by the MIT 
bag, while the binding was described by the self-consistent couplings of the confined light 
quarks (u, d) (not s nor heavier quarks) to the scalar-a and vector-o; meson fields generated 
by the confined light quarks in the other nucleons. The self-consistent response of the 
bound light quarks to the mean a field leads to a novel saturation mechanism for nuclear 
matter, with the enhancement of the lower components of the valence Dirac light quark 
wave functions. The direct interaction between the light quarks and the scalar a field is a 
key ingredient of the model, it induces a nucleon scalar polarizability 24], |25[ and generates 
a nonlinear scalar potential (effective nucleon mass), or equivalently a density- dependent 
(cr-field dependent) a-nucleon coupling. The model has opened tremendous opportunities 
for studies of the structure of finite nuclei and of hadron properties in a nuclear medium 

2ll. 



(nuclei) with a model based on the underlying quark degrees of freedom 

In QMC the Dirac equations for the quarks and antiquarks in nuclear matter, inside the 
bags of D and D* mesons, (g = u or d, and c) neglecting the Coulomb force in nuclear 
matter, are given by: 



ij-d x -(m q -V«) 1 : 1 »(V« + -V« 



H ■ d x - (m q -V*)Tl u [V« - -V« 




0. 



0. 



(16) 



(17) 

[«7 ' 9 X - m c ] ip c (x) (or ip 5 {x)) = 0. (18) 

The (constant) mean-field potentials for a light quarks in nuclear matter are defined by 
V q = g%a, VJ = g^uj and V q = g q p b, with g%., g q and g q the corresponding quark-meson 
coupling constants. 

The eigenenergies for the quarks in the D and D* mesons in units of l/R* D D * are given 
by, 




fi* ± R* DtD , [ V* + ~V* 



2 V p / 






(19) 

(20) 
(21) 



Then, the .D and D* meson masses in a nuclear medium m* D D *, are calculated by 



771 



D,D* 



2^ — 757- — + t^r dd ,b, 



]=q,q,c,c 



D,D* 



dm 



D,D* 



OR 



D,D* 



! B,B*- R D,D 



where fi* == ft£ = - [x 2 q + (R 
(R 



D,D 



*m 



1211/2 



(q = u,d), with m* = m q -g q a, tt* c 



fit 



(22) 
(23) 

[x 2 c + 



D,D 



*m. 



) 2 ] x / 2 , and Xg :C being the bag eigenfrequencies. B (=(170.0 MeV) 4 ) is the bag 
constant, n q (n q ) and n c (n 5 ) are the lowest mode quark (antiquark) numbers for the quark 
flavors q and c in the D and D* mesons, respectively, and the Zd,d* parameterize the sum 
of the center-of-mass and gluon fluctuation effects and are assumed to be independent of 
density. We choose the values (m q ,m c ) = (5, 1300) MeV for the current quark masses, and 
Rn = 0.8 fm for the bag radius of the nucleon in free space. The quark-meson coupling 
constants, g%, g^ and g q , are adjusted to fit the nuclear saturation energy and density of 
symmetric nuclear matter, and the bulk symmetry energy 21| . Exactly the same coupling 



8 



constants, g%, g^ and g q , are used for the light quarks in the D and D* mesons and baryons 
as in the nucleon. 

Because of baryon number conservation, no vector potential should contribute to the 
loop integrals. Then, the vector potentials for the D and D* mesons should be the same 
in considering the case of the DD* mixed loop to cancel out. However, for the K + meson 
case, g%, associated with the vector potential had to be scaled 1.96 times to reproduce 
an empirically extracted repulsive potential of about 25 MeV at normal nuclear matter 
density 27J. The reason is that i^-mesons may be regarded as pseudo-Goldstone bosons, 
and they are therefore difficult to describe by naive quark models as is also true for pions. 
For this reason, in earlier work we explored the possibility of also scaling the g^ strength 
by a factor 1.96 for the .D-mesons p, LZ|. In the present case, this possibility is excluded by 
baryon number conservation. As a result, the vector potential does not contribute to the 
final results. Thus, we may focus on the (scalar) effective masses of D and D* mesons. The 
QMC predictions for the in-medium effective masses of these mesons are shown in Fig. [2] 
as a function of nuclear matter density. The net reductions in the masses of the D and D* 
are near.y the same as a function of density, as dictated by ft. lu M qm rl uuraUr co mttmg 
rule [26]. 
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FIG. 2. D and D* meson (scalar) effective masses as a function of baryon density. 



IV. NUMERICAL RESULTS 



Amongst the main ingredients of the present calculation are the phenomenological form 
factors needed to regularize the self-energy loop integrals in Eq. (j7|). Following previous 
experience with a similar calculation for the p self-energy 28J, we use a dipole form for the 
vertex form factors 

:2\ 



u D ,D*{q. 



K,D* + m % 



A AD*+ 4 <D*(q) 



so that the -F](q 2 ) in Eq. ([7|) are given by 



Fdd 


(q 2 ) = 


^(q 2 ), 




Fdd 


•(q 2 ) = 


= u D (q 2 ) u D > 


(q 2 


F D * D *(q 2 ) 


= u 2 D *(q 2 ), 





(24) 

(25) 
(26) 
(27) 

where A^ and A^. are cutoff masses. Obviously the main uncertainty here is the value of 
these cutoff masses. In a simple-minded picture of the vertices the cutoff masses are related 
to the extension of the overlap region of J/ty and D mesons at the vertices and therefore 
should depend upon the sizes of the wave functions of these mesons. One can have a rough 
estimate of A^> and Ap* by using a quark model calculation of the form factors. Using a 3 Pq 
model for quark-pair creation [29] and Gaussian wave functions for the mesons, the vertex 



form factor can be written as 



4] 



Uqm{<1 



-q 2 /4(/%+2/32) 



(28) 



where (3d and (3^ are the Gaussian size parameters of the D and J/ty wave functions. 
Demanding that the it(q 2 ) of Eq. (|2~4|) and uqm (q 2 ) have the r.m.s. radii (r 2 ) 1 / 2 , with 



(r 2 ) 



-6 



dlnu(q 2 



dq 2 



i 2 =0 



one obtains 



A 2 = 32 (P 2 D + 2/3. 



Am 



D- 



(29) 



(30) 



Using m D = 1867.2 MeV and for the /3's the values used in Ref. J4|, /3 D = 310 MeV and 
(3^ = 520 MeV, one obtains Ap = 2537 MeV. Admittedly this is a somewhat rough estimate 
and it is made solely to obtain an order of magnitude estimate, since we do not expect that 
Gaussian form factors should be very accurate at high q 2 . In view of this and to gauge 



10 



uncertainties of our results, we allow the value of Ad vary in the range 1000 MeV < Ad < 
3000 MeV. Moreover, for simplicity we use An = Ad*- 

Using rriD* = 2008.6 MeV for the average of the vacuum masses of the _D*'s, there remain 
to be fixed the bare J/\l/ mass mS, and the coupling constants. The bare mass is fixed by 
fitting the physical mass mj/^ = 3096.9 MeV using Eq. (|6]). For the coupling constants 
we use g^DD = QipDD* = 9^d*d* = 7.64, which are obtained by invoking vector- meson- 



dominance and use of isospin symmetry 



30). 



We are now in a position to present the results for the in- medium mass shift Am of J/\l/, 
defined in Eq. ([5]). We calculate the in- medium self-energy using the in-medium D meson 
mass as given by the QMC model presented in Section IIHI We present results for £ = (no 
gauge coupling) and for £ = 1 (with gauge coupling). 



Ad 


m *j/i> 


DD 


DD* 


D*D* 


Am 


1000 


3081 


-3 


-2 


-11 


-16 


1500 


3079 


-3.5 


-2.5 


-12 


-18 


2000 


3077 


-4 


-3 


-13 


-20 


3000 


3072 


-6.5 


-5 


-12.5 


-24 



TABLE I. In-medium J/\£ mass th*,,^ and the individual loop contributions to the mass difference 
Am at nuclear matter density, for different values of the cutoff A^, and using the non-gauged 
Lagrangian - £ = in Eq. ([S]). All quantities are in MeV. 

Initially we present results for £ = 0. In Table |T]we present the in-medium J/\l/ mass rn*j,^ 
and the individual loop contributions to the mass difference Am at nuclear matter density 
Po, for different values of the cutoff mass Ad. First of all, one sees that the net effect of the 
in-medium mass change of the D mesons gives a negative shift for the J/\l/ mass. The total 
shift ranges 16 to 24 MeV at normal nuclear matter density. The results show in addition 
that the D*D* loop gives the largest contribution of the three. Also, this contribution is 
rather insensitive to the cutoff mass values used in the form factors. A negative self-energy 
means that the nuclear mean field provides attraction to J/^f. The important question is 
of course whether such an attraction is enough to bind J/ty to a large nucleus. A partial 



answer can be obtained as follows. One knows 



31] that for an attractive spherical well of 



radius R and depth V , the condition for the existence of a nonrelativistic s-wave bound 



11 



state of a particle of mass m is 



2fe2 



v n > 



7T 2 ti 

8mR 2 ' 



(31) 



Using for m 



m 



J/* 



and i? = 5 f m (radius of a medium-size nucleus), one obtains V > 



1 MeV. Therefore, the prospects of capturing a J/^> if produced almost at rest in a nucleus 
are quite favorable. 

Next, we assess the impact of using a gauged Lagrangian for the DD loop on m*j,^ and 
Am. The results are shown in Table HT1 The contribution of the DD loop to Am is still much 
smaller than the DD* and D*D* contributions, but of opposite sign. The net J/^> mass 
shift is still sizable, varying from 13 MeV to 18.5 MeV as the cutoff is varied from 1000 MeV 
to 3000 MeV. The small, positive value of the DD loop contribution is in agreement with 
the result of Ref. [4]. 



Ad 


m j/* 


DD 


DD* 


D*D* 


Am 


1000 


3084 


+1 


-2 


-12 


-13 


1500 


3082 


+1 


-2.5 


-12.5 


-14 


2000 


3080 


+1 


-3 


-14 


-16 


3000 


3078 


+0.5 


-5.5 


-13.5 


-18.5 



TABLE II. Same quantities as in Table [TJ but using the gauged Lagrangian - £ = 1 in Eq. (JSj) . 

In Figs. |3] - M we show the separate contributions of the DD, DD* and D*D* loops 
and their sum to the J/ty mass shift. As the cutoff mass values increase in the form 
factors, obviously each loop contribution becomes larger since the integral is divergent, but 
the increase is less pronounced for the D*D* loop. Since the D*D* loop gives the largest 
contribution, it is encouraging that this loop contribution is rather insensitive to the cutoff 
mass values used. 

V. CONCLUSIONS AND PERSPECTIVES 



We have used an effective Lagrangian approach to evaluate the D and D* loop contribu- 
tions to the mass shift of J/^> in cold nuclear matter. Effects of the medium on the D and 
D* are calculated using the QMC model, in which effective scalar and vector meson mean 
fields are coupled to the light u and d quarks in the charmed mesons. There are no free 

12 



Q 
Q 



-10 



1 

Q 

Q - 
9< 



-20 





1 1 1 1 1 


' 1 ' 1 


- 


v ^^ ' 


Trr ^: 




^ 


"- - * «T^ — — 


- 


* *s 


** ~" - -"J " * " " 


- 








DD loop 




* *-«. t 


cutoff 1000 McV 






■ ■ ■ cutoff 1500 McV 
cutoff 2000 McV 








■ — ■ cutoff 3000 McV 




1,1,1, 


1 





0.5 



1.5 



2.5 



P B /P (p =0.15fm") 



FIG. 3. Contribution from the DD loop to the difference of the in-medium and vacuum J/^f 
masses Am as a function of nuclear matter density for different values of the cutoff mass A75 . 
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FIG. 4. Contribution from the DD* loop. See also caption of Fig. [3l 

parameters in this QMC calculation since all quark-meson coupling constants and bag pa- 
rameters are fixed by fitting saturation properties of nuclear matter. The J/^f — D coupling 
constants are taken as determined from vector meson dominance and the cutoff masses are 
varied over a large range of values. The QMC predicts a 62 MeV mass drop for the D and 
D* mesons at nuclear matter density. This mass drop leads to a corresponding in-medium 
J/ty mass shift varying between —16 MeV and —24 MeV for cutoff masses within the range 
of 1000 MeV and 3000 MeV. Such a mass shift is large enough to bind a J/^f to a nucleus 
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FIG. 5. Contribution from the D*D* loop. See also caption of Fig. [3l 
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and vacuum J/^f masses Am as a function of nuclear matter density for different values of the 
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for a J/ty produced at low momentum in the rest frame of the nucleus. 

Although the conclusions of the present calculation are very promising towards the possi- 
bility of binding J/ty in a nucleus, some issues clearly require further investigation. Amongst 
the most important ones are the calculation of effective J/ty potentials for finite nuclei [32 1 
and their momentum dependence, and the inclusion of D and D* widths. Recent calcula- 
tions 33( of in-medium D and D* widths based on meson-exchange models have obtained 



14 



somewhat contradictory results and further study is required. As emphasized in Refs. [34 1 
the lack of experimental information on the free-space interaction of D mesons with nucle- 
ons is a major impediment for constraining models and the use of symmetry principles and 
exploration of the interplay between quark-gluon and baryon-meson degrees of freedom is 
essential in this respect. Still another issue is the dissociation of J/\l/ in matter by colli- 
sions with nucleons and light mesons. This subject has been studied vigorously in the last 
years using different approaches, like meson exchange [35f] and quark models |36|, QCD sum 



rules 



y 



and the NJL model 



381 ] . Finally, we stress the need for a deeper understanding of 
the role played by color van der Waals forces in the J/^> mass shift, particularly in respect 
with nucleons interacting in a nucleus. 
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